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Abstract 

Using Hamilton-Jacobi theory, we develop a formalism for solving semi- 
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Q>^ ■ classical cosmological perturbations which does not require an explicit choice 
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''si" . of time-hyper surf ace. The Hamilton-Jacobi equation for gravity interacting 

o> 

. with matter (either a scalar or dust field) is solved by making an Ansatz which 

Q . includes all terms quadratic in the spatial curvature. Gravitational radiation 

c/3 . and scalar perturbations are treated on an equal footing. Our technique en- 

compasses linear perturbation theory and it also describes some mild nonlinear 

effects. As a concrete example of the method, we compute the galaxy-galaxy 

correlation function as well as large-angle microwave background fluctuations 

for power-law inflation, and we compare with recent observations. 
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I. INTRODUCTION 



The freedom in choosing a time-hypersurface in general relativity is sometimes viewed as 
a curse because it leads to numerical problems, difficulties in quantization, etc.. However, for 
semi-classical analyses based on the Hamilton- Jacobi (HJ) equation, it is actually a blessing. 

Even in classical general relativity, selection of a useful time foliation is often a difficult 
task. In solving Einstein's equations using the ADM form (see, e.g., ref. [|T|), arbitrary 
choices must be made for the lapse N and the shift A^' functions; these fields reflect our 
liberty in choosing the time hypersurface as well as the spatial coordinates. Hamilton- Jacobi 
theory provides an elegant method of bypassing these very difficult decisions. 

It is remarkable that the Hamilton- Jacobi equation for general relativity refers neither 
to the lapse nor to the shift functions. As a result, a solution for the generating functional 
S is valid for all choices of the temporal and spatial coordinates. The HJ equation is 
the natural starting point for a hypersurface and gauge-invariant analysis — it yields a 
covariant formulation. It is analogous to the Tomonaga-Schwinger equation that was applied 
successfully to quantum electrodynamics [0. However, solutions to the HJ equation are 
difficult to obtain because one must solve for the entire ensemble of evolving universes that 
is described by superspace. 

In a series of papers 0] - we have developed a systematic method of solving the HJ 
equation by using an expansion in spatial gradients. Using the ffist few terms, one can derive 
the nonlinear Zel'dovich approximation ^ and its higher order corrections which describe 
the formation of pancake structures in a dust-dominated Universe 0-0. Causality is 
maintained in the relativistic theory, and the final expressions are actually simpler than 
those obtained from the Newtonian theory (see, e.g., Moutarde et al [T^ and Buchert and 
Ehlers |n[). 

Moreover, Parry, Salopek and Stewart derived a recursion relation which enables one 
to compute the higher order terms for the generating functional from the previous orders. 
Two useful techniques were employed: a conformal transformation of the 3-metric as well as 
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a line-integral in superspace. Deforming the contour of integration corresponds to choosing 
an alternative time-integration parameter. If the endpoints of integration were fixed, all such 
contours gave identical results provided the theory was invariant under reparametrizations 
of the spatial coordinates. Time-reparametrization invariance was closely related to spatial 
gauge-invariance . 

However, in the semi-classical problems of interest to inflationary cosmology, a finite 
number of terms is insufficient. In this paper, we will consider all terms which are quadratic 
in the curvature. We will effectively sum an infinite number of terms in the spatial gradient 
expansion. Our procedure is analogous to that developed by Barvinsky and Vilkovisky [12 



in a very different context: the one- loop effective action for gravity interacting with matter. 
In this way, we recover the results of linear perturbation theory considered during the 1982 



Nuffield workshop [|T^ - jT^, as well as those of Mukhanov, Feldman and Brandenberger 



T8| . Our formalism also describes some mildly nonlinear effects which we hope to apply to 



stochastic inflation |[T^ - Although there have been some interesting proposals, the 



choice of time-hypersurface in stochastic inflation still requires further clariflcation . For 
example, Linde et al have pointed out that eternal inflation may appear differently on 
various time-hypersurface slices. A covariant formulation would be advantageous. 

In Sec. II, we set forth the HJ equation and the momentum constraint equation. We 
consider two case of physical interest: (1) a scalar fleld in an inflationary universe and (2) 
dust in a matter-dominated epoch, with or without a cosmological constant. One can factor 
out the effects of the long-wavelength background by using a conformal transformation of the 
3-metric. We review the long-wavelength formalism using the concept of field-space diagrams 
which provide a simple illustration of the concept of hypersurface transformation. These 
diagrams are similar in spirit to Minkowski diagrams that proved useful in understanding 
special relativity. 

In Sec. Ill, we suggest an Ansatz which is second order in the spatial curvature. Arbitrary 
coefficients appear that are functions of the matter fleld and the spatial Laplacian operator. 
Substitution into the HJ equation lead to two linear differential equations, which describe the 



evolution of the scalar modes as well as the tensor modes of the 3-metric. For an inflationary 
cosmology, we choose initial conditions that are consistent with the Bunch-Davies vacuum 



There is growing interest in the gravitational waves produced during inflation because 
it was pointed out that they could provide a large part of the signal detected by the 
DMR (Differential Microwave Radiometer) experiment on the Cosmic Background Explorer 
(COBE) satellite and yet be consistent with structure formation. In our approach, 



density perturbations and gravitational waves are treated on an equal footing. In Sec. V, 
we compute large angle microwave background fluctuations for various inflationary models 
and we compare with COBE's two-year data set p8|, p9 |. 



After the ground work had been laid by Dirac (see, e.g., ref. 0]), the Hamilton- Jacobi 
equation for general relativity was first written down by Peres in 1962. For flat space- 



time, Kuchar [^] solved for the semi-classical wave-functional describing the ground state. 
By making explicit gauge choices, Halliwell and Hawking gave approximate results 
for the wavefunctional during the inflationary epoch. Salopek, Bond and Bardeen 



quantized the system using a Heisenberg formulation, and they computed the power spectra 
numerically for numerous models utilizing either one or two scalar fields. By expanding the 
action to second order in perturbations, Mukhanov et al gave an alternative prescription for 
quantizing this system in arbitrary gauges. Their final results were elegant, but because they 
perturbed the lapse and shift functions, their method was somewhat tedious. Our method 
removes this unattractive feature and it generalizes the method of Kuchar and Halliwell & 
Hawking. 

(Units are chosen so that c = SttG = Svr/mp = % = 1. The sign conventions of Misner, 
Thorne and Wheeler M will be adopted throughout.) 
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II. THE HAMILTON-JACOBI EQUATION FOR GENERAL RELATIVITY 



For a single scalar field interacting with gravity, the HJ equation and the momentum 
constraint equation are: 

ss ss 

= nix) = 7-^/^^^ [2luixhA^) - 7.-(x)7.K^)] + 
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(2.1a) 



= n,{x) = -2 U,-^] + ^^7,,^, + . (2.1b) 

V hkj{x) J J d'jkiix) d(f){x) 

In the ADM formalism, the line element is written as 

ds'^ = g^^dx^dx" = (-AT^ + f^NiNj'^ df + 2Nidt dx' + -fijdx' dx^ , (2.2) 

where and Ni are the lapse and shift functions respectively, and 'jij is the 3-metric. In 
eq.(|2.1aD, R denotes the Ricci scalar of the 3-metric. The object of chief importance is 



the generating functional S = S['~fij{x),(j){x)]. For each universe with field configuration 
['~fij{x),(f){x)], it assigns a number which can be complex. The generating functional is the 
'phase' of the wavefunctional in the semi-classical approximation: 

^ ~ e*-^ . (2.3) 

For the applications that we are considering, the prefactor before the exponential is not 
very important, although it has interesting consequences for quantum cosmology ||3^. The 
probability functional, 

p = (2.4) 



is just the square of the wavefunctional (see, e.g., ref. [§§]). It is the focus of attention 
in cosmology during inflation, and even during the matter-dominated era. The Hamilton- 
Jacobi equation (p.la| ) and the momentum constraint ( p.lb| ) follow, respectively, from the Gq 



and Einstein equations with the canonical momenta replaced by functional derivatives 
of S: 



^'\x) = -^, and vr^a;) = -^. (2.5) 
6-fij{x) d(f){x) 



Eq. ( |2.1a] ) is the relativistic generalization of the Newton-Poisson relation, whereas eq. ( |2.1b|) 



demands that the generating functional be invariant under an arbitrary change of spatial 
coordinates (see, e.g., Misner et al [|1|, p. 1185). If the generating functional is real, the 
evolution of the 3-metric for one particular universe is given by 

whereas the evolution equation for the scalar field is 

(0-iVV,)/iV = 7-^/'^. (2.6b) 

Here | denotes a covariant derivative with respect to the 3-metric 7^-,. The lapse and shift 
function appear neither in the HJ equation ( p.la|) nor in the momentum constraint (|2.1bD . 



Hence in HJ theory, all gauge-dependent quantities appear only in the evolution equations, 
( [^.6a| ) and ( |2.6bD , for the metric and scalar field. 



HJ methods can also be applied fruitfully to systems of perfect fluids Q], [^. For 



example the Hamiltonian and momentum constraints for collisionless, pressureless dust pS 
are given by: 

= n{x) = 7"^^^T — r^T — l^luixhjkix) - -fij{x)-fki{x)] 

6S 1 



+ A+r^X^ ^ - ^I'^'R + Vo , (2.7a) 

^^'^^^ + ^^^'^M + ^TT^'- (2.7b) 
^lkj{x)Jj S-fki{x) 6x{x) 

A cosmological constant term denoted by Vq has also been included. The dust field x 
describes, for example, cold-dark-matter particles. Its evolution equation is 



(x - N\,) /N = /r + ^^, (2.8) 



whereas that for the metric is given by eq.( p^.6aD . The 4- velocity of the dust is the 
4-gradient of the potential 



(2.9) 



where g'^'^ is the inverse of the 4-metric. 



A. Review of Long- Wavelength Solution 

1. Scalar Field and Gravity 

For fields where the wavelength is long compared to the Hubble radius, one may safely 
neglect second order spatial gradients (terms within square brackets) in the HJ equation 
( p.lal ) for gravity interacting with a scalar field. The resulting equation has the trivial 



solution 



5W[7.,(a;),0(x)] = -2 J d'x^'/' H [<p{x)] . (2.10) 

provided that the Hubble function H = H{(j)), a function of a single variable, satisfies the 
separated Hamilton- Jacobi equation of order zero [Q: 



Eq.( |2.10D is also a solution of the momentum contraint (|2.1b|) because the volume element 



^3^,^1/2 -g invariant under reparametrizations of the spatial coordinates. H = H((f)) corre- 
sponds to the Hubble parameter in the long-wavelength limit. We will examine in detail the 
special case of infiation with an exponential potential 

V(0) = Vo exp (-^0) , (2-12a) 

where p is a constant that describes the steepness of the potential. An exponential potential 
arises naturally in the induced gravity model (see, e.g., ref. [0) as well as in extended 



infiation |^. For this case, one can find the exact general solution of eq.( |2.11| ). In 
particular, the separated HJ equation of order zero has the attractor solution 

Vo 



1/2 



3(l-l/(3p)) 



exp = . (single scalar field) (2.12b) 



It describes power-law inflation where the scale factors evolves as a{t) oc t^, with t being a 
synchronous time variable. (In order that the scalar field convert its energy into radiation and 
matter at the end of inflation, the scalar field should have a minimum in its potential. Hence, 
we interpret eq.( |2.12a]) as describing the asymptotic branch of the potential as ^ — oo; 
we will assume that microwave background fluctuations as well as fluctuations for galaxy 
formation are generated on this branch.) The separated Hamilton- Jacobi equation has been 
widely applied in the reconstruction of the inflaton potential from cosmological observations 

El- 



2. Dust Field and Gravity 

For the case of a dust field, the long-wavelength theory is found by dropping R and 
in the HJ equation ( |2.7a| ). One can then attempt a solution analogous to eq. ( P^.IOQ except 



now H = H{x) is a function of x satisfying 



2_ 2dH ^ zt2 _ ^0 



The general solution is 



H{x) = Hq cotanh 



ix - x) 



(2.14a) 



2 

where x is a homogeneous constant 0. If the vacuum energy density Vq is negligible, then 
one recovers the Hubble parameter for a matter dominated universe in the limit that the 
cosmological constant vanishes, Hq ^ 0: 

mx)-^y (2.14b) 

The scale factor, 'j^^^, defined to be the sixth root of the determinant of the 3-metric, 

7^/^ = ^(-)(^)"'', (2.15) 

can be found by taking the derivative with respect to the parameter x 0- Here d{x) is an 
arbitrary function of the spatial coordinates. For dust with cosmological constant, we see 
that 
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2/3 



7^/6 cx d{x) \smh (^^ix - X))\ , (2.16a) 
whereas for pure dust, we recover the result 

7^/^cxci(x)(x-x)'/'- (2.16b) 
Eq.(|l6D is plotted in Fig. (2). 



B. Field-Space Diagrams 

The consequences of the long-wavelength approximation can be illuminated using field- 
space diagrams where one plots the metric variable versus the matter field. 

1. Scalar Field and Gravity 

For simplicity, we will assume that the shift function vanishes. The evolution equations 
( p.6a| ), ( p.6b| ) in the long- wavelength limit then become 



7,,/iV = 2if(0)7,j, (2.17a) 
flFf 

0/iV = -2— . (2.17b) 
These can be simplified by defining the field, 

a{t,x) = [In j{t,x)]/6, (2.18a) 

and then letting 

7.,(t,x) =e2"(*'") h,{x), (2.18b) 

where kij{x) is independent of time (this is the most interesting case) with det(A;) = l. Hence 
one obtains 

a/N = H{(f)) , (2.19a) 



0/iV = -2— . (2.19b) 

We can integrate these equations by utilizing our freedom in choosing time. If we choose 
t = to be the time variable, the lapse is defined through eq. (|2.19b| ) 

dH 

1/N = -2— . (2.20) 



Eq.( p.l9aD then becomes 

da 1 H{(p) 



2 2^ 



dH ' 



(2.21) 



which involves only a and the independent variable (p. It may be integrated immediately 
leading to 

1 f'f' . ..Hi 



a{(j), x) = ao{x) - - d(f) (2.22) 

where ao{x) is an arbitrary function of x. For the example of inflation with an exponential 
potential, eq.( p.l2b|) , we obtain the trivial solution 

a{(j),x) = ao{x) + ^(f). (2.23) 

In general, the metric variable a is inhomogeneous on a surface of uniform scalar field. For 
a single scalar field, C/3 is defined to be the metric fluctuation on a uniform slice: 

C/3 = Aa{(j)) = a{(f), X2) - a(0, Xi) = ao{x2) - «o(a^i) ; (2.24) 

it is the difference of a between two spatial points X2 and xi on a time hypersurface of 
uniform |jT^, 0], 0]. It is independent of time. For a single scalar field in the long- 



wavelength limit, this is true in general (see eq. (|2.22|) ), and not just for the example of an 
exponential potential. 

We plot the solutions eq. ( p.23| ) in the field-space diagram Fig.(l). Each curve represents 
the evolution of the fields (0, a) for a given spatial point. One may invert eq. ( |2.23| ) to obtain 
as a function of a 
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2 2 

(f)(a,x) = —\ -ao(x) + \ -a . (2.25) 
VP VP 

We interpret this inversion as choosing a time-hypersurface where a is uniform. In this 
situation, the scalar field is inhomogeneous. The fiuctuation in the scalar field between two 
spatial points X2 and Xi is given by 

[2 [2 
A(f){a) = <^(a, X2) - <^(a, Xi) = ~J ~ ["0(3^2) - ao{xi)] = ~\\~ ^C({(p) ■ (2.26) 

It is related to the metric fiuctuation through the negative of the slope of the a versus (j) 
trajectories. The transformation from a surface of uniform to one of uniform a is simply 
visualized in the field-space diagram. 

2. Dust Field and Gravity 

A similar long-wavelength analysis can be repeated for the dust field. The evolution 
equations are 

a/N = H{x) , (2.27a) 

X/N=l. (2.27b) 

If X is taken to be the time hypersurface, then = 1, and eq.( p.27a| ) may be integrated 
using the same method that was employed for a scalar field. More simply, one can just apply 



eg. (2. 16a) 



2 

a(x, a;) = ao{x) H — Insinh 
3 



SHoX 



(2.28) 



2 

We have set x = 0, which may always be arranged by shifting x- We plot a as a function 
of X ill the field-space diagram. Fig. (2). Qualitatively it is the same as Fig.(l), except that 
here the trajectories are not straight lines. 

HJ formalism is particularly useful for problems involving several fields. For example, 
in the long- wavelength limit, we have shown how to solve exactly the case of two fiuids, a 
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dust field and a field describing blackbody radiation [0. This solution describes adiabatic 
as well as isothermal perturbations. In this paper, we will consider only one matter field 
at a time. (Using a Taylor series expansion in synchronous time, Comer et al |^2[ have 
investigated various perfect fiuids.) 

These very simple considerations illustrate the role of time in general relativity. The HJ 
formalism appears to indicate that the most useful choices for the time hypersurface will be 
one of the matter fields or the metric variable a, or some combination of the two. A hy- 
persurface transformation amounts to slicing a field-space diagram in a particular direction. 
Such simple behavior is also manifested in the higher order solutions to the HJ equation for 
general relativity that we will consider in the next section. 

III. QUADRATIC CURVATURE APPROXIMATION FOR GRAVITY PLUS 

SCALAR FIELD 

A. Factoring Out the Long- Wavelength Background 

Before we attempt to solve the HJ equation ( p.la| ), we subtract out the long- wavelength 
background from the generating functional: 

S = + ^ , = -2 J £x^'I^H{ct)) , (3.1) 



where H satisfies the separated Hamilton- Jacobi equation (p.ll|) . The functional for fiuctu- 
ations, JF, now satisfies 



Although superficially similar, this step differs in principle from the usual analysis of pertur- 
bations on an homogeneous background , |^ , ]18| . Here a long- wavelength background 
is allowed, which is closely related to what is done in stochastic infiation. The first term 
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S^^^ is explicitly invariant under reparametrizations of the spatial coordinates — gauge- 
invariance is manifestly maintained. Moreover, no explicit choice of the time parameter has 
been made. 

The first line of eq.( |3.2|) may simplified if one introduces a change of variables, [jij, (p) — *• 



U 



where the conformal factor Vt = VL{u) is defined through 

d\nVL _ ^dHd\nVL_^ 

du dcj) d(j) 

Functional derivatives with respect to the fields transform according to 
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(3.3a) 



(3.3b) 



(3.4) 



6<p 2^\6u-^^ '-5h, 
In order to simplify the notation, we will henceforth suppress the symbols |„ and \f.. which 
denote the variables that are held constant during differentiation. Utilizing the conformal 
3-metric fij instead the original 3-metric is analogous to using comoving coordinates 
rather than physical coordinates in cosmological systems. At long-wavelengths, a surface 
of uniform u corresponds to comoving, synchronous gauge because = 1 in eq. (|2.19b| ). 
Even if (considered as a function of H) oscillates, u is monotonic. However, when one 
considers short- wavelength terms associated with the functional T for fluctuations, a surface 
of uniform u no longer represents a synchronous gauge. 
The HJ equation reduces to 



5£_ 

6u 



2H{u)U, 



-\ 2 



5/, 



6u 



+ n-'(v) f-V2 ^^^^ i2f f f f 1 + ^''(^)^''^' 

^^"^ ^"^^ jf-.Ji;^^'^''^^'-^-^'^^^ 8(f)^.) 

where the functional S^"^^ is given by 

S^^\fij{x),u{x)]= J d^xf^^[j{u)R + k{u)u'\i] . (3.6a) 

From now on, a semi-colon will denote a covariant derivative with respect to the conformal 
3-metric, e.g., = f^-^u-j. In addition, R is the Ricci scalar of fij. The w-dependent 
coefficients j and k, 



(3 
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j(u) = / du' +F, k{u) = H{u)n{u) , (3.6b) 

Jo 2 

where F is an arbitrary constant, were first derived in refs. [Q, ||^ in order that the spatial 
gradient terms appearing in eq.(|3.2|), 



1 , 1 , ^5(2) 

(3.7) 

may be expressed as a functional derivative with respect to u{x) holding fij{x) fixed. The 
momentum constraint maintains the same form as before but it is now expressed in terms 
of the new variables {fij,u). 

n.(x) = -2(u^] +^f,,^^+^-fu, = 0. (3.8) 

B. Integral form of HJ equation 

It will prove useful to work with an integral form of the HJ equation. Before proceeding, 
we pause to consider a simple illustration from potential theory. 

1. Potential Theory 

The fundamental problem in potential theory is: given a force field g^{uk) which is a 
function of n variables u^, what is the potential $ = ^{uk) (if it exists) whose gradient 
returns the force field: 

^=9'{uk) ? (3.9) 

OUi 

Not all force fields are derivable from a potential. Provided that the force field satisfies the 
integrability relation, 

" d d 



duj dui 



duj dui 



$, (3.10) 



(i.e., it is curl-free), one may find a solution which is conveniently expressed using a line- 
integral 
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Huk)= / Y.dvj9'{vi)- (3.11) 

If the two endpoints are fixed, all contours return the same answer. In practice, we will 
employ the simplest contour that one can imagine: a line connecting the origin to the 
observation point Uk- Using s, < s < 1, to parameterize the contour, we define the 
variable vi = vi{s,ui), 

vi = sui with dvi = ds ui (3.12) 

so the line-integral may be rewritten as 

" /■! r/s 

Hu,) = J2 -v.g^ivk) . (3.13) 
j=l -^0 s 

Similarly, in solving for the generating functional, we will employ a line-integral in su- 
perspace. The integrability condition for the HJ equation follows from the Poisson bracket 
of the constraints (Moncrief and Teitelboim |^ ) 

{n{x'),n{x''')} = [Y'{x')nj{x') + Y^{x'')n,{x'')] ^^s\x'-x''). (3.i4) 

In fact, each contour of the line-integral corresponds to a particular time-hypersurface choice. 
Provided that the generating functional is invariant under reparametrizations of the spatial 
coordinates, (e.g., Tij vanishes in eq.( p.l4D ), different time-hypersurface choices will lead 
to the same generating functional. Hypersurface invariance is closely related to gauge- 
invariance. 



2. Line-Integral in Superspace 



Following Parry, Salopek and Stewart 0, an integral form of the HJ eq.( |3.5| ) may be 
constructed using a line-integral in superspace: 

ds „ , . . „ _]^i2 bT 



T\UM),u{x)\^ \ d^x \ ^vn-'{v)r 

J Jo s 



d X 



1 dsvQ-^{v)f-^/^ 



s 
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[^fafjk — fijfki] + 



1 2 



5f, 



*j J 



S^'^Mx),uix)] (3.15) 



In analogy with the example from potential theory, we replace the index j with the spatial 
coordinate x and the finite sum J2j with the integral / d^x. The integrating parameter is 
again denoted by s, and v{x) = su{x) represents a straight line in the superspace of the 
scalar field; however, at each spatial point, fij{x) is held fixed in the line-integral. As a result, 
one may safely add an arbitrary functional of fij to the right-hand side. No approximations 
have been made thus far; eq.( p.l5|) is exact. 



C. ANSATZ 

For a spatial gradient expansion of S we have already shown how to set up a recursion 
relation which gives higher order terms from the previous orders. Explicit solutions 
were given which were accurate to fourth order in spatial gradients. However, we now wish 
to investigate higher order terms. Finding an explicit and exact expression which is valid 
to all orders appears to be extremely difficult. Instead we will examine an infinite subset 
consisting of all terms quadratic in the Ricci curvature, but containing any number of spatial 
gradients. Halliwell and Hawking used a similar approximation, although our expansion 



will manifestly maintain hypersurface and gauge-invariance. We hence make an Ansatz of 
the form, 

jr = 5(2)_^g^ ^ith (3.16) 



Q = f d^xf/^ \r S{u, D^) R + R'^ f{u, D^) Rij - Ir f{u, D^) R 
J V 8 . 



(3.17) 



where S{u, D"^) and T(m, D'^) are differential operators which are also functions of u. is 
the Laplacian operator with respect to the conformal 3- metric, e.g., 

D^R = D''D,, R = R',i = f-^'^ [f'^PRj) ^ . (3.18) 

We refer to Q as the quadratic functional. We interpret the operator 5* for scalar perturba- 
tions to be a Taylor series of the form 
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S{u,D')=Y.SMiD'r, (3.19) 

n 

and similarly for the operator T describing tensor fluctuations. Note that T is sandwiched 
between two Ricci tensors, Rij, as well as two Ricci scalars R. The full Riemann tensor 
Rijki does not appear in the Ansatz because for three spatial dimensions because it may be 



written in terms of the Ricci tensor (see, e.g., [0). Although the flrst and third terms in 



eq. (|3.17|) may be combined into a single one, the present form simplifies the final evolution 



equations for S and T (see eqs.( |3.27a| ,b)). 



1. Order of Perturbation 

We now define the rules which determine the order in perturbation of various terms. By 
first order, we refer to terms such as R or D^u, or D^R or D'^u, which vanish if the fields 
are homogeneous; they may contain any number of spatial derivatives. Quadratic terms are 
a product of two linear terms. Some examples are u-iU'^, D^uD^R. It should be clear how 
to determine other cases. 

Other quadratic terms such as {D'^uY and RD'^u, could be included in the Ansatz (p.l7|) 



as well. However for the case of a single scalar field as well as that for a single dust field, it 
may be shown that they actually vanish. This is one of advantages of utilizing the conformal 
transformation ( |3.3a| ,b). (This simplification does not occur for multiple matter fields, and 



analogous terms would have to be included.) 

All that is necessary now is to compute some functional derivatives, and substitute into 
the integrated HJ equation (|3.15| ). It is useful to note that for a small variation of the 



conformal 3-metric 6fij the corresponding change in the Ricci tensor is 

^^ij ^ 2'^'^' [^fli;jk + ^flj;ik — ^ fij;lk " Sflk;ij] ■ (3.20) 

In the integral form of the HJ equation, integration by parts is permitted which simplies the 
analysis considerably. (However, in its differential form ( |2.1a| ) one cannot simply discard 
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total spatial derivatives.) In addition, all cubic terms are neglected. For example from 
eqs.( p.64 b) we find that 



(3.21) 



5u 5fij 

where a term proportional to u-iU'^ has been dropped, otherwise if ( ^.21| ) where squared, this 
term's contribution to the integral HJ equation ( p.l5|) would be cubic and higher which we 
are not considering here. Other terms that appear are 



d-^x I — v^l -^{v) f 



1 ds 

s dfij dfki 

^ n/ 



[^filfjk — fijfkl] 

3 



/ d'xj^ ^/V2 2vn-%v)j\v) [Rm.,, - ^B?) . 



(3.22) 



The order of some operators may be changed if the commutator is an undesired higher order 
term, e.g.. 



(3.23) 



Hence, effectively commutes with any function of u. For a similar reason, one may 
permute the order of spatial differentiation with impunity, D-iD-j ~ D-jD-i, because the 
additional Riemann tensor term would be of higher order. As a result, we find 



d X / — V 



-1/2 



s 



(f(«))' 



ov 



d^xf/'^R / dv\ r 

^0 \8n%v){^{v)y 



6iH\v) S^D^ + 16H{v) 



^2,„,^ c27^4 , ,«^,.,^/^^v-; H{v) j{v)] SD' + 



-H{v)j{v) 



>R, (3.24) 



and that 



/ d'x f ^ ^^-V-^/^lf If [2fuf,k - Mm] 
J Jo s dfij dfki 



ij <JJkl 



J d^xf^uR^£dvn- 



2f + 4j TD'^ + 2 T'^D^ 



Ri 



—R / dvn-^ 

8 io 



2f + Aj TD' + 2 T^D^ 



(3.25) 



Collecting terms, the integral form of the HJ equation ( p.l5| ) gives: 



0= [£xf/^R S+ r dv ^ ^ {^^-Hiv)]iv) + 9>H{v)SD^ 



T+ dv 
Jo 







Rjj — —R 

^ 8 


f + / 


JO 



R + 



r\. (3.26) 



Each coefficient must vanish separately, leading to a pair of uncoupled integral equati 



ions. 



Taking the derivative with respect to u for each, one obtains the differential equations 

(3.27a) 



0- — + 

du 



Since one may add an arbitrary functional of fij to the integral HJ eq. ( P.15| ), these first order, 
ordinary differential equations are equivalent to the pair of integral equations arising from 
eq. (p.26| ). These nonlinear differential equations are of the Riccati type. They represent 
tremendous simplification over the original HJ equation ( p.lal ). 



a 



D. SOLVING RICCATI EQUATIONS 



As is well known, the Riccati equations, (|3.27a|) and (|3.27b|) , may be reduced to linear 
ordinary differential equations. To this aim, we define the Riccati transformation 



through 



w = w{u,D'^), y = y{u,D'^ 



8H^D^ w du 



Idy ~ ^ 



(3.28) 



2/^4 ydu 



(3.29) 
(3.30) 



which lead to the desired result: 



d'^w { d 



In 



H dd) 



— Vr'^{u)D'^w , (scalar perturbations) (3.31) 

ou 



= + 3H(u)— Vr'^{u)D'^y . (tensor perturbations) 

ou'^ ou 



(3.32) 
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Equation (|3.31|) describes the evolution of scalar perturbations (density perturbations) of the 
metric, whereas eq.( |3.32| ) describes tensor perturbations (gravitational waves). We empha- 
size that no time choice has been made in deriving these basic equations, although the scalar 
field parametrizes the evolution because u = u{(f)) is the independent variable. These are 
the fundamental equations of cosmological perturbations that we will utilize in this paper. 
If we set, 



1 H 



(3.33a) 



we obtain the equation for scalar perturbations that was derived by Mukhanov et al p8|] 

= ^ + miuV^ + [m\^^ - Q'^D^)z , (3.33b) 



where the effective mass rrieff is given by 



m 



eff 



d'^V dH 
h 2 

ddP' du 



3 - 



1 dH 

+ 2 



m du 



du^ 



H- 



du 



(3.33c) 



In deriving these equations, it is useful to note several relationships between the scalar field 
and the new variable u which follow from eqs.( p.34 b) and ( p.ll| ): 

dH 



du 

dU 



du^ 



.3H^ - — 
du dd) ' 



du^ I du du dctP' du^ I du 



(3.34a) 
(3.34b) 
(3.34c) 



Eq.( |3.34b| ) is the well-known evolution equation for a long- wavelength scalar field. It is 
identical to the evolution equation in synchronous gauge (A^ = 1) describing a homogeneous 
Friedman universe. 



1. Exact Solution 



Many exact solutions have been found for power-law inflation (see, e.g., eq.( |2.12bD ). 
Abbott and Wise have shown that the evolution equation (|3.32|) for tensor perturbations 
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can be solved exactly. In addition, Lyth and E. Stewart have shown that the same 
solution may be applied to the scalar perturbation eg. ( |3.33b"D . For example, if we use H{(f)) 
defined in eq.( |2.12b| ), we find that 

1 dH 1 , , 

(3.35) 



H d(j) y/2p 

is a constant, and the equation for scalar perturbations eq.( |3.31| ) becomes 

= — — H u w . (scalar perturbations) 

ow^ u ou 



(3.36) 



Here we have also used the fact that for power-law infiation, the Hubble parameter and Vl 
are given by eqs.( p.34 b) 



P 



H{u) = — , and ^l{u) 



u 



(3.37) 



Eq. ( p.36| ) is identical to that describing a massless scalar field in a Friedman universe. This 
is quite surprising since one would have naively expected that the mass term d'^V/dip'^ in 
eq.( |3.33(^ ) would lead to damping of the fiuctuations with wavelengths larger than the Hubble 
radius but evidently the additional terms cancel this effect. In fact, the equation ( p. 36] ) for 



scalar perturbations is identical to the tensor evolution eq.( 3.32 ). The solution for both of 
these equations can be expressed in terms of a Hankel function of the first kind: 

(3p - 1) 



'-w(u,D^ 
V 



7]-u(l-3p) 



with u 



2{p-l) 
(3.38) 



The normalization of this solution is irrelevant since a logarithmic derivative enters in the 
Riccati transformation ( ^.29] ), ( |3.30D . However, in order to agree with the conventions of 
Birrell and Davies we have assumed that 



dy dy* ' 



^%u){y*^-y 



(3.39) 



du du ^ 

Note that the argument for the Hankel function is negative, and it increases to zero as u 
varies from to oo. For an infiationary epoch where p > 1, it describes a positive frequency 
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solution as M — s> 0, e.g., when the wavelength of a particular mode is far within the Hubble 
radius, 

(-D^Y^^^ ( r-=- r du \ 

y = z^ ^ . ^ exp i^-D^ 7^] (u^O). (3.40) 



As a result, eq. (|3.38|) corresponds to the ground state wavefunctional at wavelengths shorter 



than the Hubble radius (i.e., the Bunch-Davies vacuum |2^; see also ref. ||3^). However, 



when the wavelength exceeds the Hubble radius, the state is no longer in the ground state. 
(There is also an analogous solution to eg. ( |3.36| ) which involves a Hankel function of the 
second kind but it describes a negative frequency solution at short wavelengths — it does 
not describe a system which is initially in the ground state.) 

IV. QUADRATIC CURVATURE APPROXIMATION FOR GRAVITY PLUS 

DUST FIELD 

The quadratic curvature approximation used in Sec. HI may be applied to any field that 
is derived from an action principle. We derive the scalar and tensor equations corresponding 
to a dust field. Unfortunately one cannot define a ground state for the dust field as was the 
case for a scalar field. Hence the conditions at the beginning of the dust-dominated era were 
generated at earlier times, such as during the scalar field dominated epoch of inflation. 

We follow the treatment given for a scalar field. In analogy to eq.( ^lD , we express 
the generating functional S = S[yij{x),x{x)] for dust and gravity as the sum of a long- 
wavelength part S^'^^ and a fluctuation part JF. Here, however, the Hubble function H = 
H{x) is given by eq.( p.l3| ). The HJ equation ( p. 7a]) for dust then becomes 



Since we cannot solve this equation exactly, we will retain only those terms which are at 
most quadratic in perturbation: 
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We will also employ a conformal transformation of the 3-metric —>■ fij 



dH 

dx 



(4.2) 



fij = ^ ^ix)7ij 



(4.3) 



where fl = ^{x) is defined through 



dx 



H{X). 



(4.4) 



(For example, if the vacuum energy density Vq vanishes in eq.( p.l3| ), the Hubble function H 
and the conformal factor can be written as 



H = —, n = x^'\ (withx = 0) 
3x 



(4.5) 



Hence, the HJ equation for dust becomes 



OJij Ofkl 



[2fiifjk — fijfki] 



(4.6) 



fij 



where the functional S^"^^ is given by 



5(2)[/,^.(x),x(x)] = J d'xf/' [j{x)R + k{x)x''X-4 ; (4.7a) 

once again, a covariant derivative with respect to fij is denoted using a semi-colon. The 
X-dependent coefficients j and k are f^, [§, 

rx nix') 



Jix) = r dx' + F, k{x) = H{x) Q{x) , 
JO 2 



(4.7b) 



where F is a constant, so that 



(4.8) 



may be expressed as a functional derivative with respect to x{^) holding fij fixed. Using a 
line-integral in superspace, we can construct the integral form of this equation analogous to 
eq.(CT: 
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Jo s dfij dfki 



(4.9) 



where v{x) = sx{x) represents a straight hne in superspace. This equation has a similar 
but simpler form than the case of a scalar field. We make the analogous quadratic Ansatz 
as in eq. ( |3.17| ), and we find that the evolution equations for the scalar and tensor operators, 
S = Six, and f = fix, are 

dS 



0= 
0= 



dx 
df 



^ (jix)+fD^'" 



dx ^%X) 

The solution to the scalar perturbation equation ([4.10a|) is trivial: 



(4.10a) 
(4.10b) 



S = SiD'^) (scalar perturbations) 



(4.11) 



is an arbitrary functional of the conformal Laplacian operator. It is independent of the dust 
field X- Since dust never oscillates like a scalar field, one cannot identify a ground state 



for the dust field. The tensor perturbation equation (|4.10b|) may be solved using the same 
technique that was employed for a scalar field. We first define y, 

1 dy 



T 



which leads to a linear equation. 



2D^ ydx 



-^-D-'jix) 



(4.12) 



Q'^y (}y ' 

= 7— r + SHix)^^ ^ '^ix)D'^y , (tensor perturbations) 

dx dx 

describing the evolution of the graviton in a universe with dust. 

If the cosmological constant vanishes, we obtain the exact solution: 



(4.13) 



y = Ax-^'^Ji 



Bx-"'J- 



(4.14) 



The choice of the coefficients A and B will be given in the next section. In fact, the forms 
for S and T at the beginning of the matter-dominated epoch are determined by a preceding 
period of inflation with a scalar field. 
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V. LARGE ANGLE MICROWAVE BACKGROUND FLUCTUATIONS AND 

GALAXY CORRELATIONS 



We first show how to compute large angle microwave background anisotropics arising 
from both the scalar and tensor fluctuations of inflation. In principle this is not difficult 
since it essentially amounts to expanding in spherical harmonics. We first interpret the 
semi-classical wavefunctional that was computed in the previous sections. 

By maintaining gauge and hypersurface invariance in the dynamical analysis, we have 
not made any extraneous assumption which would typically complicate the results. Con- 
sequently, the final equations for a scalar field, ( p.31| ), (|3.32|) and for a dust field, ( [4.10aD , 
( [4.13| ) are of a very simple form. However, many measurement processes actually choose a 
specific time-hypersurface. A simple example arises in special relativity: when measuring 
the lifetime of a particle, one typically chooses the rest frame of that particle. In an appli- 
cation of high interest to cosmology, Sachs and Wolfe derived the large angle microwave 
background anisotropy by utilizing comoving, synchronous gauge (uniform x slice). 



A. Interpretation of Semi-Classical Wavefunctional 

In earlier papers 0] - 0, it was shown that the gradient expansion could be used to 
compute nonlinear effects in cosmology. In particular, we computed higher order corrections 
to the Zel'dovich approximation. However in computing microwave background anisotropics, 
it is sufficient to consider a linear approximation; we will now consider only a small deviation 
hij of the conformal 3-metric fij from fiat space 5ij , where the probability functional reaches 
its maximum: 

fij{u, x) = 6ij + hij{u, x) with hij{u, x) = J2 77^-^3 <^'^'^ ^aiu, k) Elj\k) . (5.1) 

We have expressed hij as a sum of plane waves with comoving wavenumber k. Furthermore, 
we also expanded it in a complete basis of 6 symmetric, polarization matrices E\^\k). For 
example, if the comoving wavenumber is aligned with the z-axis, we choose them to be: 
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E^^^ = Diag[l,-1,0] , (tensor) (5.2a) 



= = 1 , (all other components vanish) (tensor) (5.2b) 

4^ = Diag[l,l,0] , (scalar) (5.2c) 

Eff =Diag[0,0,v^], (gauge) (5.2d) 

e[^^ = E^f = 1 = £'23'' = -^32'* • (s-ll other components vanish) (gauge) (5.2e) 

The first two are traceless and divergenceless; they correspond to tensor perturbations. The 
third describes scalar perturbations, whereas the remainder are gauge modes. The matrices 
have been normalized so that 

= 26'''' . (5.3) 
In addition, one must respect the reality condition: 

K{k) = Ki-k) . (5.4) 



To linear order, the Ricci tensor for the conformal metric may be computed using eq.( p.20|) 
with fij = 6ij and 6fij = hij to give 

Rij = ^ {hii/ + hij/ - hij/ - h^ij^ . (5.5) 

The probability functional arising from inflation is then given by the square of the wave- 
functional, eq.(|2r 



m'hA^). m] = exp - - / \\(3,{k)\' + \(32{k)\' + \(3s{k)\'] , (5.6a) 



2 J (27r) 



where 



h,ik) = ^ A(fc), hik) = ^ h{k) , h^{k) = ^ /?3 

^2k^Ti{u, -F) ^2k^Ti{u, -P) ^16k^Si{u, -k^) 

(5.6b) 

Here has been replaced by —k'^. In Fourier space, S and T are complex numbers which 
we expand into real and imaginary parts: 
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S{u,-P) = SR{u,-k^)+iSj{u,-k^), and f = Tr{u, -k^) + iTj{u, -P) , (5.7) 

with 



1 (^Y 1 1 1 



which follows from the Riccati transformation equations (|3.29|) , ( p.3C ) and the normalization 



conditions (|3.39|) ; (we have assumed that u is uniform (see below) and j(M) is real). Because 



the wavefunctional is invariant under reparameterizations of the spatial coordinates, the 



gauge modes hi{k)^ h^{k) and hQ{k) are absent in eq. (|5.6a|) . Hence they are unrestricted, 
and they may assume arbitrary values consistent with the reality condition ( |5.4| ). The three 
/3's are Gaussian random fields which satisfy the following: 

< I3a{k) f3;{k) >= (27r)3 6'{k - k') 6ab , Mk) = f3:{-k) , a, 6 = 1, 2, 3. (5.9) 



The polarization matrices have been chosen so that the probability functional ( |5.6a| ) is 
diagonal in hi{k), h2{k) and h^{k). (Using a canonical transformation in conjunction with 
HJ theory, Langlois |^ has derived a reduced phase space Hamiltonian; our approach 
differs from his in that we perform a phase space reduction after finding a solution to the 
HJ equation.) 

1. Long-wavelength Fields from Power-law Inflation 

Long-wavelength fields are measurable and their evolution was discussed in Sec. H.A. 
Recall that we defined C/3 to be the fluctuation in a = (ln7)/6 on a comoving slice, 
eq. (|2.24|) . By choosing u = u{(f)) to be uniform in eq. (|5.6a|) , we see that for small deviations 
from flat space that zeta is related to h^lk) through 

C(x) = -y J^s^""'h3{k). (5.10) 
Since the gauge modes are arbitrary, we have chosen = h-^/y^ in order that 

h{u, k) E§\k) + Kiu, k) E\f{k) = h^iu, k)6,j (5.11) 
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is proportional to the identity matrix. It has become conventional to define the power 
spectrum through 

fc3 



27r2 



d X e 



-ik-x 



<ax)m>. 



(5.12) 



which leads to 



27r2 



9\w\ 



(5.13) 



Loosely speaking, [k^lwl"^ / {2tt^)Y^^ can be interpreted as the metric fluctuation, Aa{(f)) = 
C/3, on a uniform slice (see eq.( |2.24| )), whereas [fc^|z|2/(27r2)]^/2 jg ^j-^g fluctuation in 
the scalar fleld, A0(a), on a uniform a (uniform curvature) slice (see eg. ( p.26|) ). They 
are related through the hypersurface transformation eq.( p.33a| ); see Sec. II. B and Fig.(l). 
(In fact, Hwang has used uniform curvature slices to provide an elegant derivation of 
the scalar perturbation equation (|3.33b|) .) In DeSitter space where the scale factor varies 
exponentially in synchronous time, Q{u) = e^", where if is a constant, it is useful to note 
that that at long-wavelengths 



.1/2 



(27r2 



(27r2 



1/2 



H 
2^ 



(5.14) 



is given exactly by the Hawking temperature [^, |3J|. If inflation is not exactly DeSitter, 
this is approximately true in which case it is useful to interpret H = H{k) as the Hubble 
parameter at the time that the comoving scale crossed the Hubble radius during inflation, 
k/{Hn) ~ 1, and 
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(5.15) 



For power-law inflation, we will compute this quantity exactly. 

For the tensor modes, we deflne hi{x) and h2{x) as well as their corresponding power 
spectra, e.g., 

d^k 



hi{x) 



(5.16) 
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with 

V,Ak) = 1^ / d'x e-^'-^^ < h{x)h{0) >=^2\y\\ (5.17) 

Using the exact solution eq.( p.38|) in the long- wavelength limit, u — >■ oo, we find for power- 
law infiation that y and z are independent of time u: 

-(3p-l) 

The long-wavelength power spectra for power-law inflation are power-laws which are related 
to each other through the steepness parameter p, eq. ( p.l2a| ) : 

\ rts — 1 

— j , (long - wavelength) (5.19a) 
4 ( kV~^ 

VhM = VhM = ^ ^c(^o) ( ^ j , (long - wavelength) (5.19b) 

where, in the case of power-law inflation, the spectral indices, ns,nt, for scalar and tensor 
fluctuations actually coincide 

n, = nt = l- . (5.20) 

ip-l) 

V(^{ko) is the value of the power spectrum in zeta at some flducial wavenumber scale fco 
which we will choose later. (It has been pointed out that inflation with a cosine potential 



can also yield a power-law fluctuation spectrum for scalar perturbations |51 . 



2. Heating of the Universe 

At the end of inflation, the scalar fleld typically rolls to the minimum of its potential 
where it oscillates and converts its energy into a thermal bath consisting of radiation and 
matter. Because the coupling of the scalar fleld to matter is not well understood, several 
possibilities can arise yielding either a low or high value of T^ax, which we deflne to be the 
maximum temperature reached immediately after inflation. (Many models utilizing super- 
symmetry give values, Tmax ~ lO^GeV, whereas the Variable Planck Mass model ||34 
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is rather special in giving a very high result, Tmax ~ lO^^GeV.) Fortunately, the amplitude 
of the metric at large wavelengths is indifferent to the uncertainty in T^ax- (However, the 
present physical length of a comoving scale is indeed sensitive to the value of the maximum 
temperature). Numerical calculations for a simple phenomenological model describing cou- 
pling of a single scalar field to radiation and matter demonstrate that the amplitude of the 
metric fluctuation ( on a. comoving slice remains constant ||3^. Essentially a result of mo- 



mentum conservation, the rest frame of the single scalar field (uniform 0) slice is coincident 
with the rest frame of radiation which is identical to the rest frame of the matter (uniform x 
slice) at wavelengths larger than the Hubble radius. Hence the fluctuations for structure for- 
mation arising from inflation are typically adiahatic. (If there are two scalar fields which are 
important during inflation, one may also produce isocurvature perturbations which we will 
not consider here; see, e.g., Sasaki and Yokoyama [Q.) Moreover, the function y describing 
the tensor fluctuations in eqs.( |3.30| ) and ( ^4.12|) is continuous during the heating process. 



If we are only interested in large angle microwave background fluctuations produced in 
the cold-dark-matter model, we may thus equate the probability functionals on comoving 
time slices before and after heating of the Universe: 

VYiij\u = Uheat] = 'P[lij\X = Xheat = O] . (5.21) 

Hence for long-wavelength fields in the radiation and matter dominated eras, the power 
spectra for zeta as well as that for tensor perturbations are identical to those arising from 
inflation, eqs. (|5.19a]) and ( |5.19b|) . In other words, the scalar and tensor operators, S and T 



are continuous on a comoving time slice. 



B. Application of Sachs- Wolfe Formula 

It remains useful to continue employing comoving slices of uniform x since the phase 
transition where radiation becomes uncoupled from matter occurs on a slice where the tem- 
perature is uniform, T = 4000i^, which for adiabatic fluctuations, coincides with a uniform 
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X slice at large wavelengths (see, e.g., ref. [0). Hence by observing the phase transition, 
one has effectively chosen a very special time-hypersurface. The Sachs- Wolfe formula 
yields the large angle temperature anisotropy from a line integral over the perturbation in 
the 3- metric hij computed in uniform x gauge [comoving, synchronous gauge): 

Ar(x)/T= -ieV£rfx|^[x,x(x)] (5.22) 

The line integral traces the path x{x) of a photon path from the surface of last scattering to 
the present epoch; e* is a unit vector giving the direction of the photon's velocity. For angles 
of interest to COBE (a > 7°), we are concerned with those comoving scales that reenter the 
Hubble radius during the matter-dominated era. 

1. Scalar Perturbations 

In comoving synchronous gauge, the scalar part of the metric evolves according to 

h^\r,x) = '^ax)6,, + ^r2 0,(x) , (with N = 1,N' = 0). (5.23) 

(This may be derived most simply by using ref. P], or by applying Sec. IV.) Here we have 
defined conformal time r = t{x) to be 

r.pf.^^(Xy\ aM^Uf', (5.24) 
Jo a{x) Ho \xoJ \XoJ ^Hq 

where Hq is the present value of the Hubble parameter which we assume to be 50 km/s /Mpc 
which is consistent with measurements of the Sunyaev-Zel'dovich effect by Birkinshaw et al 
|5^ (see also Lasenby |5^). In order to agree with our present units of measurement, we 
have normalized the scale factor a(x) oc ^ix) to unity at the present epoch xo- Integrating 
eq.( |5.22| ) twice by parts and retaining only those boundary terms which are important for 
a > 2*^, one finds that the large-angle temperature anisotropy, 

Ar(x)/T = -C(x)/15, (5.25) 
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in a flat, matter-dominated Universe is proportional to the value of ({x) on the surface of 
last scattering which is a sphere of radius R = \x\ = 11, 700Mpc. In this way, temperature 
fluctuations are a probe of scalar fluctuations from inflation. 

Since we are concerned with temperature anisotropics measured on the celestial sphere, it 
is natural to employ a spherical harmonic expansion. A plane wave can be decomposed into 
orthogonal spherical harmonic functions, Y^m, and spherical Bessel functions, je, through 
(see, e.g., ref. |^ ) 



^ik-- = 4^ ^ ^ i^j^^kr)Y,l{n^) F,™(fifc) . (5.26) 

£=0m=-£ 

If 6 denotes the angle between k and x, then the addition theorem relates the spherical 
harmonics with the Legendre polynomial Pe{cos6): 

Piicos6) = T.yLm Yimi^k) . (5.27) 

Hence, the plane wave expansion 

ax) = / -— e^'- — Vdk) (3,{k) , (5.28) 



(27r)3 \^ P 

(which follows from the expression for the scalar power spectrum eg. ( |5.12|) of a Gaussian 



random fleld) implies that the scalar contribution to the angular correlation function Cs{q) 
can be expressed as a sum over the Legendre polynomials P^(cosa) ||57[| : 

oo 

Cs{a) =< AT{x)AT{x') >,= ^ < AT| >, Pi{cosa) , (5.29) 

e=o 

< AT^ (^^j (2^ + 1) ^Vdk)jKkR) (5.30a) 

^ ^2 (2^ + 1) r(£+(n.-l)/2) r((9-n,)/2) 

5 T{i+{5-ns)/2) r((3 + r2,)/2)' ^ ' ^ 

A-T^V,{koKkoR) -2 p(2-n./2) r((9 - n.)/2) ^^'^^"^ 
where a is the angle between the two points, x and x', on the surface of last scattering; = 
2.736 ± 0.0177^^ is the mean background temperature [0]. Vc_{k{)) was deflned in eq.( p.l9a| ), 
and we will assume that the flducial wavenumber scale is k^ = 10~^Mpc. (Integrals of 
various combinations of Bessel functions may be found in Gradshteyn and Ryzhik |]59|). 
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2. Tensor Perturbations 



The derivation of microwave anisotropies from tensor perturbations is similar in principle 
to the scalar case. However it technically more complicated because the angular correlation 
function obtains contributions from points within the surface of last scattering. 

For primordial gravitational waves described by eq.( ^.19bD , the metric for the tensor 
modes evolves according to 

fcg' (X. = / e'-- (^n. (*)) ^ t ^S?' (*) ■ (5-31) 

during the matter-dominated era. By applying a continuity argument after the exit of 
inflation to the radiation and matter dominated eras, we have determined the coefficients 
A and i? = in eq. (|4.14| ); we have also chosen to rewrite eg. ( |4.14| ) in terms of a spherical 



Bessel function jiikr) of order 1. Using the Sachs- Wolfe formula eg. ( |5.22| ), we find that the 
temperature anisotropy is 

(5.32) 

and the resulting two-point correlation function is 



a=l 



where 

describes the paths of the photons, with e* = xY|a;|, e'* = x''^/\x'\. The values of conformal 
time at the present epoch and at the time of decoupling, (1 + a)^^ ~ 1300, are, respectively, 



T2 = 2/Hq and ti = 2/ (i/o vlSOO). This expression may be simplified by noting that 

2 

^ E^f (k) E^'^\k) = 6ip6jq + 6ig6jp- 6ij6pg , where 6ij = 6ij - kikj/k"^ , (5.35) 



a=l 
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so that 

2 

e'e^e'^e"' Eff{k) E^^^{k) = 2 (cosa - cos^cos^')^ - (l - cos^^) (l - cos^^') . (5.36) 



a=l 



Here 6 denotes the angle between k and x, and analogously for 6', and once again, a is the 
angle between x and x'. The various factors of cos^ (and cos6'') may be removed by several 
applications of the identity, 

(^ + 1) 



cos^^ Pz-fcos^ 



-P^+i(cos^) + 



-P£_i(cos6 



(5.37) 



(2£+l) ' (2£+l)' 

in conjunction with the plane wave decomposition into Legendre polynomials, eqs.( ^.26"D and 
( ^.27| ). The angular integrations in eq. (|5.33|) may be performed, and once again the angular 
correlation function can written in a series of Legendre polynomials. 



9, 



< AT'f >t= ^T^ (i - !)£{£ + 2){2£ + 1) 

dww V,, [^^) [J^ dv 
cos(w(l — f)) sm{w{l — v)) sm{w{l — v)) 



ji+2iwv) 



w^(i — vy 

ji{wv) 



+ 



(5.38) 

(5.39) 
(5.40) 
(5.41) 

(5.42) 
lO-^Mpc"^ 



_(2£ + l)(2£ + 3) (2£-l)(2£ + 3) (2£+ 1)(2£- 1)_ 
which was derived by Abbott and Wise (see also Starobinsky 0|). kr, 
corresponds to the Hubble radius at matter-radiation equality. For angular scales that are 
measured by COBE {a > 7^), the precise value of kmax is not that important, since the 
spherical Bessel functions je{ujv) provide a cutoff. 

Since the scalar and tensor contributions are independent Gaussian random fields, they 
add in quadrature 



< ATf >=< ATf >, + < ATf >t 



(5.43) 



In Fig. (3), we have computed the relative contribution of the tensor component < AT/ >t 
I (< AT/ >s + < AT/ >() for various values of the spectral index n = Ug = rit = 1— 2/(p— 1). 
For smaller values of ra^, it increases quite dramatically. 
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C. Recent Observations 



The COBE DMR team has recently analyzed their 2- year data set pM]. At the 68% 



confidence level, they find that the spectral index |29] for scalar perturbations is = 
1.10 ± 0.32, which is consistent with the simplest models of infiation models which yield 
Ug < 1. They also determine that the root-mean-square temperature anisotropy with dipole 
removed is askyiW^) = 30.5 ± 2.7 [iK at the same level of confidence. The latter quantity, 

oo 

f^sfcy(10°) = E < > exp [-/(/ + 1)/13.52] , a,fcj,(10°) = 30.5 ± 2.7^K , 

1=2 

is computed using the sum of the scalar and tensor fiuctuations eq. (|5.43|) . It determines the 
arbitrary normalization factor appearing in long- wavelength power spectra, eqs. (|5.19a] ,b). 
The exponential factor corresponds to a Gaussian window function with full width at half 
maximum of 10*^. 

However, COBE by itself cannot discriminate between tensor and scalar fiuctuations. 
One needs an additional experiment to measure the scalar perturbations. Two proposals 
have been suggested using either: (1) galaxy clustering data ||26| or (2) intermediate mi- 
crowave background experiments 1° < a < 2° |61[]. We shall discuss only the first proposal 



here since there are large variations in the intermediate angle observations ^2[ - ||65|| . 

In Fig. (4), we have computed the power-spectra for zeta that arises from the vari- 
ous power-law infiation models, assuming that they account for COBE's measurement of 
crsfcy(10°). In the limit that ri^ ^ 1 (p — oo) gravitational waves do not contribute to 
COBE's signal, and the power spectrum for zeta is the fiat, Zel'dovich spectrum. As Ug 
decreases, gravitational waves are significant, leaving a smaller contribution for the scalar 
perturbation to crsfcy(10°). Hence at the fiducial wavenumber ko = 10~^Mpc~^ (scales probed 
by COBE), V(^{ko) decreases as decreases. Moreover, the slope of the power spectrum 
becomes more negative as decreases. 

In comoving synchronous gauge (uniform x), p7^/^ is independent of time and the linear 
density perturbation at early times is 
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^2 

6{t, x) = {p-p) /p = -— Cj{x) , (early times) (5.44) 

which may be derived from the expression for the metric eq.( p.23| ). In Fourier space this 
yields, 

6{T,k)='^ak)- (5.45) 

During the radiation-dominated era, density perturbations oscillate and they damp because 
of the Hubble expansion. This effect is described by the transfer function T{k), so that in 
the matter-dominated era the linear density perturbation is given by 

5(r,fc) = ^T(fc)C(fc). (5.46) 

This requires an assumption for the dark matter, and we have adopted the cold-dark-matter 
transfer function . 



In Fig. (5), we show power-spectra for the density perturbation. 



Vs{t, k) = ^J d'x e-*^' < <5(r, x)6{t, 0) > , (5.47) 
arising from power-law inflation. The bold line depicts the power-spectrum 

Vt''\k) = - sin (^] r(2 - 7) {kroT (5.48) 



TT V 2 

for the observed galaxy-galaxy correlation function ^gg{r) = (r/rg)"'^ where tq = 10 Mpc, 
7 = 1.8, and F is the gamma function. The corresponding biasing parameter bp is found by 
computing the mass fluctuation on a scale of 16 Mpc: 

/AM \ ^ 

< (r = 16Mpc) >=l/bl. (5.49) 



V M 

For Us = 1,0.95,0.9,0.85,0.8,0.7,0.5, we compute the biasing parameter to be bp = 
0.82,1.06,1.34,1.65,2.0,2.88,5.6. In order to be consistent with the biased galaxy forma- 
tion, we insist that bp < 2 and 0.8 < Ug < 1 otherwise there are not enough fluctuations to 
seed galaxies. As a result, for power-law inflation no more than 50% of COBE's signal can 
arise from gravitational waves. 
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Previously, it had been suggested that power-law inflation with = 0.5 pH] could ac- 



count for the APM (Automatic Plate machine) survey ||6^ which demonstrated more power 
than predicted by the standard cold-dark-matter model. However, Salopek pointed out 
that the production of gravitational radiation, which was neglected in the previous calcula- 
tion, could be quite signiflcant for power-law inflation. He was able to rule out this promising 
model for large scale power since it is essential that > 0.8. (For a careful discussion of 
statistical limits on the spectral index Ug using only COBE data, consult Kurki-Suonio and 
Mathews pj). 



VI. CONCLUSIONS 

Hamilton- Jacobi methods are a cornerstone of modern theoretical physics, and they 
may be profitably applied to numerous problems in cosmology. For example, they have been 
successfully employed in deriving the Zel'dovich approximation and its higher generalizations 
from general relativity [0 ^ @]- Various researchers have employed HJ methods in an attempt 
to recover the inflaton potential from cosmological observations ^Ij. Moreover, they can 
be used to construct inflationary models that yield non-Gaussian primordial fluctuations 
0]; such models could possibly resolve the problem of large scale structure |^. Here we 
have given a careful and detailed computation of the galaxy-galaxy correlation function and 
large-angle microwave background fluctuations arising from power-law inflation, which is 
the most interesting model involving gravitational radiation. We flnd that the resulting 
spectral index for scalar perturbations must satisfy Us > 0.8, otherwise the production of 
gravitational radiation is excessive, and there are not enough fluctuations to seed galaxies. 

Our analysis is greatly facihtated by the fact that a choice of the time-hypersurface is 
not required in the computation of the probability functional during the inflationary epoch. 
Field-space diagrams are useful in visualizing a hypersurface transformation. However in 
the end when one compares with observations one typically assumes a particular choice of 
gauge. For large-angle microwave background fluctuations which are computed using the 
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Sachs- Wolfe formula, comoving synchronous gauge is preferred. 

Our line integral formulation of the HJ equation ( p. 151 ) goes a long way in illuminating 
the role of time in semi-classical general relativity. Different time-hypersurface choices cor- 
respond to different choices of contours in superspace. Provided spatial gauge invariance is 
maintained, they all yield the same result for the generating functional. 

A complete quantum formulation of the gravitational field is still lacking. String theory 
is a possible candidate, and its applications to cosmology are currently being investigated 
|70| , |7l| . Our aim is more modest in that we have restricted ourselves to the semi-classical 



theory of Einstein gravity which is nonetheless adequate in describing various quantum 
gravitational phenomena including graviton fluctuations beginning initially in the ground 
state ||2^. We hence follow in spirit the historical development of the theory of atomic 



spectra. Before the development of the quantum theory in 1926, the semi-classical theory 
of Bohr and Sommerfeld provided a useful although imperfect description of various atoms. 
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VII. FIGURE CAPTIONS 



Fig.(l): In the long- wavelength hmit, the evolution of the scalar field and the metric are 
shown for inflation with an exponential potential. For each spatial point corresponds a 
trajectory (which is a straight line). Slicing this diagram in a particular direction would 
represent making a time-hypersurface choice. A0(q;) denotes the scalar fleld fluctuation on 
a time-hypersurface of uniform a = (ln7)/6. Aa{4>) refers to the metric fluctuation on a 
time-hypersurface of uniform 4>. These perturbations are related to each other through the 
slope of the a — (f) trajectories. 

Fig. (2): In the long- wavelength hmit, the evolution of the dust fleld and the metric are 
shown when a cosmological constant is present. Aq;(x) refers to the metric fluctuation on a 
time-hypersurface of uniform x (comoving, synchronous gauge). Ax(q;) denotes fluctuation 
of the dust fleld on a time-hypersurface of uniform a. The hypcrsurface transformation 
relating the two is more complicated than in Fig.(l) because here the trajectories are curved. 
Fig. (3): The gravity wave contribution to large angle AT/T can be dominant for power- 
law inflation which employs an exponential potential, V{(f)) = Vq cxp[— As a function 
of the spherical harmonic i, the relative contribution of the gravity waves < AT| >t to 
the total contribution < AT^ >s -|- < AT/ is plotted for various potential parameters 
Us = 1 — 2/ (p—l). If one normahzes to COBE, then Ug > 0.8 is required to give fluctuations 
large enough to produce galaxies. 

Fig. (4): Primordial scalar perturbations of the metric are described by the function (. The 
fluctuation spectra for zeta are shown for various choices of the the spectral index Us arising 
from power-law inflation. They have been normalized using COBE's 2-yr data set. 
Fig. (5): For the present epoch, the power spectra for the hnear density perturbation S 
in comoving synchronous gauge are shown. The dark line depicts the observed two-point 
correlation function describing galaxy clustering. If there is no biasing, = 0.9 gives a good 
flt to the observed data near k — 10~^ Mpc~^. In order that there be enough fluctuations to 
seed galaxies, one requires that the biasing parameter bp be less than 2 which implies that 
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0.8 < < 1. As a result, for power-law inflation, at most 50% of COBE's signal may arise 
from gravitational waves. 
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FIGURES 




Fig.(l) 



FIELD-SPACE DIAGRAM: DUST FIELD + GRAVITY 




Fig. (2) 



INFLATION WITH AN EXPONENTIAL POTENTIAL 
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POWER SPECTRA FOR ZETA 
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